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Abstract 



We report a detailed study of the electronic and structural properties of the 
39K superconductor MgB2 and of several related systems of the same fam- 
ily, namely Mgo.5Alo.5B2, BeB2, CaSi2 and CaBeo.5Sii.5. Our calculations, 
which include zone-center phonon frequencies and transport properties, are 
performed within the local density approximation to the density functional 
theory, using the full-potential linearized augmented plane wave (FLAPW) 
and the norm-conserving pseudopotential methods. Our results indicate es- 
sentially three-dimensional properties for these compounds; however, strongly 
two-dimensional cj-bonding bands contribute significantly at the Fermi level. 
Similarities and differences between MgB2 and BeB2 (whose superconduct- 
ing properties have not been yet investigated) are analyzed in detail. Our 
calculations for Mgo.5Alo.5B2 show that metal substitution cannot be fully 
described in a rigid band model. CaSi2 is studied as a function of pressure, 
and Be substitution in the Si planes leads to a stable compound similar in 
many aspects to diborides. 

PACS numbers: 74.25.jb, 74.70.-b, 74.10.+V, 71.20.-b 
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I. INTRODUCTION 



The recent discovery of superconductivity at ^ 39K in the simple intermetalhc com- 
pound MgB2 is particularly interesting for many reasons. In the first instance, this critical 
temperature is by far the highest if we exclude oxides and Cgo based materials. Furthermore, 
B isotope effect indicates that MgB2 may be a BCS phonon-mediated superconductor, 
with Tc above the commonly accepted limits for phonon-assisted superconductivity. An- 
other reason of interest is that the AlB2-type structure of MgB2 (wich can be viewed as an 
intercalated graphite structure with full occupation of interstitial sites centered in hexagonal 
prisms made of B atoms) is shared by a large class of compounds (more than one hundred) , 
where the B site can be occupied by Si or Ge, but also by Ni, Ga, Cu, Ag, Au, Zn and 
Cd, and Mg can be substituted by divalent or trivalent sp metals, transition elements or 
rare earths. Such a variety is therefore particularly suggestive of a systematic study. To 
confirm this interest, a recent study |0] showed that superconductivity is supressed by ad- 
dition of more than about 30% of Al atoms. Pressure effects also lead to a depression of 
sup er conduct i vity. 

While claims of very high temperature superconductivity in diborides were never 
confirmed, superconductivity at Tc ~ 14K was observed in disilicides. Recent experimental 
work has in fact pointed out that pressure can be effectively used to tune both structural and 
physical properties of the metal silicon compound CaSi2 . In particular, Bordet and co- 
workers 1^ identified a new structural phase transition for CaSi2 from trigonal to hexagonal 
AlB2-type structure, with superconductivity at Tc = 14 K, one of the highest transition 
temperatures known for a silicon based compound. The occurrence of this structural phase 
transition implies that Si shifts from sp^ to sp"^ hybridization, and experiment points to a 
clear enhancement of when going from an sp^ to an sp"^ environment |^ . The temperature 
dependence of the Hall coefficient observed in CaSi2 indicates that this compound is a 
compensated metal for wich electrons and holes contribute in an equivalent way to the 
conductivity in qualitative agreement with electronic structure calculations. 

Early studies @-10| of the electronic structure of diborides by the OPW and tight-binding 
methods have shown similarity with graphite bands. A more systematic study of AlB2-type 
compounds [|T2| , pr3| by the full-potential LAPW method showed in simple metal diborides the 
presence of interlayer states, similar to those previously found in graphite [^1|,|T3|, in addition 
to graphite-like a and tt bands. Transition metal diborides, on the other hand, show a 
more complex band structure due to the metal d states. A very recent electronic structure 
calculation [|I6| on MgB2 shows that metallicity is due to B states. Phonon frequencies 
and electron-phonon couplings (at the F point) lead these authors to suggest phonon-based 
superconductivity. Further, new calculations [|r^ on related compounds {LiB2 and graphitic- 
like B) addressed the peculiar contributions to conductivity coming from the interplay of 
2D and 3D Fermi surface features. An alternative mechanism has been proposed by Hirsch 



11|, based on dressed hole carriers in the B planes, in the presence of almost filled bands. 



As for the silicides, many electronic structure studies have been carried out for the different 



phases |T8|-pO|, while the AlB2-type structure of CaSi2 was first studied by Kusakabe et al. 
pT| , before its actual synthetization. 

As previously pointed out, this family of compounds provides a wide playground for the 
search of new superconductors, and we intend to provide here a contribution, based on a de- 
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tailed study of electronic and structural properties together with phonon spectra, for several 
AlB2-type compounds and for CaSi2, in the trigonal and hexagonal structures. The ques- 
tions that we want to address are the following: (i) how do the electronic properties of MX2 
depend on the atomic species M; namely, can rigid band schemes provide the correct trend, 
or is the metal effect mostly mediated by lattice parameter changes (chemical pressure), (ii) 



Ref. shows the presence of both B vr Fermi surfaces (FS) giving rise to 3-dimensional 



states, and of small 2-dimensional B a-bonding pockets. It is of course interesting to know 
which of the FS is mostly involved with superconductivity. The trends observed in the band 
structures, together with experimental studies, will contribute to elucidate this point. In 
the same way, changes in the Fermi level density of states with chemical composition will 
possibly correlate with Tc values, (iii) Pure CaSi2 shows a trigonal instability, joined by a 
buckling of Si2 planes, which is apparently strongly reduced under high pressure, or by Be 
substituting some of the Si sites. The presence, or the proximity, of structural instabilities 
(investigated through phonon spectra), will hopefully provide hints to further experimental 
search of new materials. 

In this paper, we present calculations for MgB2, Mgo.5Alo.5B2, BeB2, CaSi2 and 
CaBexSi2-x (2; = 1); based on the local density approximation (LDA) to the density func- 
tional theory. We use both the Full Linearized Augmented Plane Wave method, and 
the norm-conserving pseudopotential method. In Sec. ^ we give some computational detail; 
in Sect. |T| we discuss the electronic properties for all the compounds considered; we finally 



draw our conclusion in Sect. [V 



II. COMPUTATIONAL AND STRUCTURAL DETAILS 

Our calculations were performed using the "all-electron" full potential linearized aug- 
mented plane waves (FLAPW) [|22| method, in the local density approximation (LDA) to 
the density functional theory. In the interstitial region we used plane waves with wave vector 
up to Kmax = 3.9 a.u. and 4.1 a.u. for disilicides and diborides, respectively (we used the 
larger Kmax = 4.4 a.u. for BeB2 because of the smaller spheres imposed by the smaller 
lattice parameters). Inside the muffin tin spheres, we used an angular momentum expansion 
up to Imax = 6 for the potential and charge density and Imax = 8 for the wave functions. The 
Brillouin zone sampling was performed using the special fc-points technique according to 



the Monkhorst-Pack scheme |£3[, and also the linear tetrahedra method up to 60 A;-points 
in the irreducible Brillouin zone. We used muffin tin radii Rmt = 2.2, L8, 1.59 and 2.1 a.u., 
for Ca, Si, B, and the remaining elements respectively (we used Rbc = 2 a.u. and Rb = 1.5 
a.u. in BeB2). 

The crystallographic structure of AIB2 (C32) is hexagonal, space group P6mmm with 
graphitic B planes, and 12-fold B-coordinated Al atoms sitting at the center of B hexagonal 
prisms. The unit cell contains one formula unit, with Al atom at the origin of the coordi- 
nates and two B atoms in the positions (di = (a/2, a/2\/2, c/2), d2 = — di). The trigonal 
distortion in CaSi2 corresponds to a buckling of Si planes leading to a P3ml space group 
symmetry ), with one internal order parameter z - being z = 0.5 the value corresponding to 
the AIB2 structure. 



Phonon frequencies were computed using the ABINIT code by means of Linear 



Response Theory P3|J2H] (LRT) approach in the framework of pseudopotentials plane wave 
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method p^. Troullier and Martins soft norm-conserving pseudopotentials were used for 
all the elements considered, Mg and Ca p semicore states were not included in the valence, 
their contribution to the exchange and correlation potential was accounted by the non- 
linear core corrections [^]. The local density approximation to the exchange-correlation 
functional is used. A plane wave cutoff of 44 Ry was found sufficient to converge structural 
properties of the compounds. Brillouin zone integration was performed using a 12 x 12 x 12 
Monkhorst-Pack mesh. To improve fc-point sampling convergence a gaussian smearing 
with an eletronic temperature of 0.04 Ha was used. Experimental lattice constant were used 
in the ground state and LRT calculations. 



III. ELECTRONIC AND STRUCTURAL PROPERTIES 

A. Diborides 

We shall start our discussion of electronic properties of diborides by presenting in Fig. |1| 
the energy bands of the superconducting compound MgB2, which will be the reference for 
all our further studies (the M — T — K lines are in the basal plane, while the L — A — H lines 
are the corresponding ones on the top plane at = tt/c). We used the experimental lattice 
parameters a = 3.083A, c = 3.52A. Our results, which are in excellent agreement with those 
of Kortus et. al. |]16[, show strong similarities with the band structures of most simple-metal 
diborides []T2|, computed by the FLAPW method. They show strongly bonded sp"^ hybrids 



laying in the horizontal hexagonal planes, forming the three lowest (a-bonding) bands, as 
the main structure in the valence region. The corresponding antibonding combinations are 
located around 6 eV above the Fermi level (Ep). 

The bonding along the vertical direction is provided by the vr bands (at ~ -3 eV and ~ 
2 at r), forming the double bell-shaped structures. Similarly to AlB2-type diborides, and 
unlike graphite, we notice a large dispersion of the vr bands along the kz direction [T — A). 
This can be understood on a tight-binding approach: the phases of the Bloch functions 
lead the B p^ orbitals on adjacent layers to be antibonding at F, and bonding at A. The 
downwards dispersion of the vr-bonding band (from —3 to —7 eV along T — A) is joined by 
an opposite upwards dispersion of an empty band, (at 2 eV at the F point) wich is a direct 
indication of the strong interaction among these states. The analysis of the wavefunctions 
for this empty band reveals that this band has an interstitial character 64% of its charge 
is located in the interstitial region), and is very similar to those found in graphite [|15 



m 



diborides [|T^], in CaGa2 []T3[, but also in CaSi2. In fact, this band is missing in atomic- 



orbitals only based tight-binding calculations |]T0[, and has been shown |jT2[ to be present 



even in hypotetical compounds completely deprived of intercalated metal atoms. A general 
description of the bonding can be given by the partial density of states (PDOS), shown in 
Fig. 1^. The B PDOS, in particular, shows bonding and antibonding structures for the s and 
p states. The Pz states give rise to a rather wide structure responsible for less than half of 
the DOS at Ep, which therefore results having predominantly B p — a character. 

In the search of similar materials having the desired superconducting properties, BeB2 
is the first natural candidate, as the band filling is expected to be pretty similar, and lighter 
Be atoms may help providing larger phonon frequencies while keeping similar electronic 
properties. The only experimental report p9| on BeB2, to our knowledge, gives the average 
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lattice constants. These parameters have been used for all the previous calculations of the 
BeB2 |^,12|- We have optimized the lattice constants, obtaining a = 2.87A, c = 2.85A, 



smaller than in MgB2 and similar to the values estimated by the average experimental 
values. In particular, we notice that while a reduces only by about 7.4%, c reduces by 
about 24% when Mg is substituted by Be, a result which can be understood pretty easily. 
On one hand, in fact, the optimization of a-bonds prevents a too drastic variation of the a 
parameter; on the other hand, c can change more easily as the vertical bonding, provided 
by the vr bands, has a large contribution from the metal orbitals. The band structure of 
BeB2, shown in Fig R|, are pretty similar to those computed by the FLAPW method |T^ 



at the experimentally estimated lattice parameters. We can remark the strong similarities 
with the MgB2 compound, which includes the presence of a-bonding hole pockets along 
r — ^4, but also a few relevant differences. First of all, the shorter lattice parameters lead 
to wider valence bands, and also to more dispersed a bands, in particular the cylindrical 
hole pocket along T — A, which might be relevant for superconductivity. Also, the different 
energy location of the metal s-free electron band (lower in MgB2), and the different shape 
of the TT-bonding band (related by the different c values) causes a different occupation of 
this band, especially along F — M. 

To further investigate these two materials, we studied the zone center phonon frequencies 
of MgB2 and BeB2, using the linear response approach, with the pseudopotential method; 
calculations have been carried out at the experimental lattice constant. Our results are 
listed in Table |. The two lowest independent frequencies correspond to a motion of the 
intercalated metal atoms relative to the rigid B networks: one vertical mode (^2u), and 
one doubly degenerate in-plane mode, (-Eim). The two higher frequencies represent the 
internal motion of the B network itself: again one vertical mode (-Big), now representing a 
trigonal distortion producing a buckling of the planes similar to that found in CaSi2, and 
one degenerate {E2g) mode giving an in-plane stretching of B-B bonds. Our calculations 
for MgB2 are in reasonable agreement with those of Kortus et al. |]TB[, apart from the 



frequency, higher in our calculation. This discrepancy may be due to the different 
computational schemes used, as well as to the use of the experimental lattice constant, 
probably differing more from the equilibrium value within pseudopotential than FLAPW 
method. The comparison between the two compounds shows a much higher frequency 
in BeB2. Surprisingly, however, two of the remaining frequencies are lower in BeB2, namely 
the mode which involves the motion of lighter Be atoms against the B network and the 
Big mode. This indicates a strong dependence of the phonon spectra of these materials on 
the electronic structure details. The lowest Ei^ mode has comparable frequencies. 

The states at the Fermi level, responsible for superconductivity, show two different or- 
bital characters: p — a bonding (coloumn-like FS around T — A) and pz, having vr-bonding 
and antibonding character on the basal and on the top {kz = tt/c) planes, respectively. 
Which one of them is the most important for superconductivity is not yet known. A very 
recent experimental report by Slusky et al. shows that Al doping destroys bulk super- 
conductivity when the Al content x is greater than ^0.3. Stimulated by these results, we 
studied the chemical substitution of Mg with Al. The substitution corresponding to a 50% 
concentration of Al (x = 0.5) leads to a reduction of the cell parameters, from a = 3.083 A, 
c = 3.521 A 13 for MgB2, to a = 3.047 A, c = 3.366 A 0] for Mgo,5Alo.5B2. These 



differences again indicate hardly compressible B-B a bonds, and interlayer distances quite 
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sensitive to the intercalated cation. The study of Slusky et al. p[ shows consistent results, 
and indicates, furthermore, that after a two-phases region for Al substitutions from 10 to 
25%, c collapses and bulk superconductivity disappears. We first studied this problem in a 
rigid band approach, bringing MgB2 to the same lattice parameters of Alo.3Mgo,7B2 [@. The 
only apparent change at that electron concentration is the complete filling of the tt bonding 
band at the M point. The a bonding bands need a larger electron addition to be completely 
filled in a rigid band model. 

To investigate further this problem, we studied the system with 50% Al concentration. 
We used in our calculations the experimental values, and in a first approach simulated the 
disorder (certainly present in the real compound) using an orthorhombic supercell, con- 
taining two formula units. We started by calculating the energy bands using the MgB2 
lattice constants, and then we used the experimental Mgo.5Alo.5B2 constants. The results 
of these calculations are shown in Fig. ^ together with the bands of MgB2, folded into the 
Brillouin zone (BZ) of the orthorhombic structure. The folding is such that the F point 
of the orthorhombic BZ corresponds to the F and M points of the hexagonal BZ, while Z 
corresponds to both the hexagonal A and L. Therefore, the T — Z line collects the T — A 
and the M — L hexagonal lines, and we notice immediately the three tt bands (two bonding 
and one anti-bonding, coming from folding). The T — X line contains the (folded) hexagonal 
F — M line. If we first look at the frozen structure calculation Fig. ^(c), we see quite similar 
band dispersions with the obvious splittings, related to the symmetry lowering induced by 
the different cations, and a different band filling. We notice, however, a lowering of the 
folded TT antibonding band, related to the different cations: the vr band is in fact pushed 
downwards by the interaction with the free-electron metal-s band, which is at lower energy 
in the compounds with Al (this band is partially filled in pure, non superconducting, AIB2 
|T^). As a consequence, the extra electrons coming from Al do not fill the cr-bonding band 



completely, leaving hole pockets around Z. As we relax the structure, the a bands increase 
their width due to the smaller a value, and the situation near Ep changes only because 
of the occupation of the free-electron band around F. These results suggest that although 
the rigid band scheme can be considered roughly correct, the dependence upon doping of 
the physical properties and of Fermi surface states relevant for superconductivity, may be 
crucially dependent on the actual substitutional atom. 

We have computed, within a rigid band scheme and using the scheme described, e.g., in 
Ref. the Fermi velocities, plasma frequencies (along the principal axes of the crystal), 
and the independent components of the Hall tensor for BeB2 and MgB2. Because of hexag- 
onal symmetry, the non- zero components will be Rxyz-, corresponding to the magnetic field 
along the z axis and transport in plane, and Rzxy = Ryzx corresponding to in-plane mag- 
netic field. The results, plotted in Fig. |^ as a function of the hole or electron-type doping, 
show relatively smooth variations over a wide range of doping, and are nearly coinciding, at 
zero doping, with the corresponding calculations by Kortus et al. |jl6|]. After comparing the 
two materials, the following remarks can be done: (i) the density of states at Ep is smaller 
in BeB2. (ii) MgB2 is nearly isotropic in terms of plasma frequencies and Fermi velocities, 
while BeB2 is expected to show important anisotropy in resistivity, and, in particular, the 
in-plane conductivity is expected to be smaller. A band-by-band decomposition explains 
these results in term of very similar (and very anisotropic) contributions in the two com- 
pounds from a states, and quite different contributions from the n bands, nearly isotropic 
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in MgB2 and favoring conductivity along the c-axis in BeB2. These results point out that 
despite structural similarity, the electronic properties of AlB2-type structures differ substan- 
tially from graphite and intercalation graphite compounds, (iii) The Hall coefficients are 
positive (hole-like) when the magnetic field is along the z axis (transport in-plane), while 
is negative for MgB2 and almost zero but negative for BeB2, when the magnetic field is 
in-plane. A detailed analysis shows again comparable contributions in the two compounds 
from a bands, while the sign of Rzxy and Ryzx comes from a balance of the vr bands. The 



Hall coefficient of MgB2 has been recently measured by Kang et al. |3^. Ru-, which will 
correspond in this case to an average of the tensor components, turns out to be positive 
and to decrease with temperature. At T = lOOK they give Ru = 4.1 x lO^^^m'^/C, and 
an extrapolation of their results at T — 0, gives Rh ~ 6.5 x lO^^^m'^/C. If we average, at 
zero doping, our tensor components, the positive contributions overcome the negative ones, 
resulting in an average value Rh ~ 2 x 10~^^m^/C. This value is thus of the correct order, 
but single crystal measurements are probably necessary to have a significant comparison, 
avoiding an average of quantities with different sign. It is worth noticing that, if the sign 
of carriers has a role on superconductivity ||1T|, BeB2 also has a positive Hall coefficient. 
Summaryzing this discussion, we see that if superconductivity would turn out to be related 
to (T-bonding electrons, the two compounds will be basically different (electron-wise) only 
because of the larger density of states in MgB2 which in turn can be ascribed to the different 
lattice parameters. 



B. Calcium Disilicide 

When prepared at ambient pressure CaSi2 has a rhombohedral structure (space group 



i?3m, a = 3.85 A, c = 30.62 A) |Q and is a semimetal |T9|, not superconducting down to 
30 mK 0]. By annealing under pressure - typical condictions are 4 GPa and 1300 — 1800 K - 
a tetragonal superconducting (a — ThSi2— type) phase appears with a 3-dimensional network 



3|| of Si, showing superconducivity with = 1.58 K. At pressures between ^ 7 GPa and 
9.5 GPa rhombohedral structure samples undergo a transition to a trigonal phase (P3ml, 
ar ~ 3.78 A, ct ~ cr/6 ~ 4.59 A) with Si at 2d positions, and z = 0.4. A second 
phase transition occurs at 16 GPa. This new phase may be described as [§] an AlB2-type 
structure with z slightly smaller than the ideal value 1/2. The transition to the AlB2-type 
structure implies a switch of the Si bonding from the tetrahedral sp^ to an almost planar 
sp"^ coordination. Kusakabe and coworkers [E^ have predicted the existence of an AlB2-type 
polymorph at high pressure. The AlB2-type structure was previously observed in rare-earth 
silicides [^, while in CaSi2 it is apparently stabilized by doping the Si planes with Be 



36| . We will discuss this point later. We report electronic structure calculations for the 
AlB2-type polymorph of CaSi2 at the experimental lattice costant, as given by Bordet et al. 
0], and we discuss the relative stability of the hexagonal versus trigonal polymorph, after a 
presentation of the electronic states. 

In Fig. ^ we show the energy bands of CaSi2 in the AIB2 structure (dashed lines), 
superimposed with those of the trigonal structure (full lines), along the main symmetry 
lines of the Brillouin zone. Since this polymorph can be obtained only under pressure, we 
do not minimize the lattice parameters, but rather use the experimental values a = 3.7077 A 
and c = 4.0277 A, and optimize the internal parameter z. The bands in the AIB2 structure 
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show strong similarities with those of simple metal diborides, and of CaGa2 [|T2|- The major 



difference is provided by Ca d states, which are mostly found in the conduction region, above 
4 eV; however, there is a p-d mixing between Si and Ca, especially concerning Si vr and Ca 
dz2 orbitals, also mixing with the free-electron states. Unlike MgB2, the Fermi level of CaSi2 
cuts the TT-antibonding bands over all the Brillouin zone. This of course implies a reduced 
contribution of tt states to the stabilization of an sp"^ environment, consistent with the sp^- 
like trigonal distortion. We will come back to this point later. The energy bands in the 
trigonal structure show that, a part from a small rigid shift, the a bands change very little 
with respect to the AIB2 structure. The tt bands, on the other hand, show anticrossings 
which in the antibonding case are located right at the Fermi level. This is most clearly 
indicated in the density of states of CaSi2, reported in Fig. |^. The dip at Ep in the trigonal 
structure provides a textbook explanation for the lattice distortion. While it is tempting to 
associate the enhancement in the AIB2 structure with the larger density of states at Ep 
within a BCS approach, other factors as phonon spectra and electron-phonon couplings will 
obviously be as important. 

In order to understand the stability of AlB2-type polymorph, we study the total energy 
variation as a function of z, relative to the ideal z = 0.5 case, for two experimental structures 
1^, and report it in Fig. ^(a). There is a good agreement between theory and experiment, 
and a non-negligible (~ 30 mRy per formula unit) stabilization energy. As we fix a and c to 
the values corresponding to the larger experimental pressure P], we see much smaller changes 
in the equilibrium z than experiments {z ~ 0.448), relative to the lower pressure values. This 
indicates that some different modification might be going on in the experimental samples. 
In order to assess the stability of this compound in the AIB2 structure, and on the light of 
our previous discussion, we vary the atomic number of cation, which now describes a virtual 
atom ranging from Ca to K (corresponding to Z = 20 + x). Our results, reported in Fig|^(b), 
show that while x = —0.3 leaves the total energy curve almost unchanged, x = —0.7 or larger 
give a nearly vanishing trigonal stabilization energy. This is consistent with our previous 
conclusion that a filling of vr-antibonding bands destabilizes the AIB2 structure. This last 
conclusion is supported by the existence, in the hexagonal structure, of CaBexSi2-x, with 
X = 0.75, with a = 3.94 A and c = 4.38 A. In fact. Be doping removes electrons from 
the Si2 planes, and brings the formal filling of the above vr-antibonding bands to the level 
found in diborides. To investigate these effects, we studied the (artificially ordered, but 
existing) CaBexSi2-x system for x = 1, not far from the experimental stoichiometry. The 
corresponding bands, shown in Fig. ^, have been calculated at the experimental lattice 
parameters ordering Be and Si atoms as, e.g., in hexagonal BN. In this way, the Brillouin 
zone is the same as in the undoped compound, and the interpretation is easier. The CaBeSi 
band structure shown are overall similar to that of the undoped compound. The symmetry 
lowering induced by substitution produces the large splitting of degeneracies at the K and 
H points. The major difference is of course related to the position of the Fermi level, which 
is almost exactly coinciding with that found in the superconducting MgB2. The reduced 
bandwidth for the a bonding bands may be associated with the smaller size of Be orbitals. 
The electron depletion from the antibonding vr-bands removes the trigonal instability, as 
demonstrated by total energy calculations. The Big mode, corresponding to the trigonal 
distortion, will therefore be close to an instability in Be-doped calcium disilicide, which may 
offer an interesting playground in the search of superconducting materials. 
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While it is difficult to say, on the basis of our electronic structure calculations only, 
whether the similarities with MgB2 band structure can lead to superconductivity, it may be 
useful to provide hints on the stability of this compound, which may help experimentalists 
in the search of new compounds. The PDOS in Fig. |10| show a strong hybridization between 
Si and Be states, in the valence region. In other words, substitutions on the Si site can be 
possible without disrupting the sp^ network if one-electron energies of the impurity atoms 
are comparable with those of Si, and light atoms should be preferred in order to keep high 
values of phonon frequencies. We may suggest B, which however corresponds to a hole 
doping smaller than Be. 

IV. CONCLUSIONS 

We have performed electronic structure calculations for MgB2, BeB2, Mgo.5Alo.5B2, CaSi2 
in the trigonal and in the AZi?2-type polymorphs at the high pressure experimental structural 
parameters, and for the CaBexSi2-x system (x = 1), experimentally stabilized for x = 0.75. 
We have calculated band structures, PDOS, transport properties and phonon frequencies, 
using the FLAPW and the norm- conserving ab-initio pseudopotential methods, within the 
local density approximation. The following conclusions can be drawn on the basis of our 
calculations: as compared to MgB2, BeB2 has a similar ffiling of a bands, but differs in 
terms of tt states near Ep. The DOS at Ep is lower, and some phonon frequencies are 
substantialy higher while others are comparable or lower. It would be interesting to see 
whether this compound is superconductor, but only one experimental report exists for BeB2. 
Al substitution cannot be fully simulated by a rigid band model, as the tt antibonding states 
appear to be quite sensitive to Al substitution, and absorb much of the added electrons. 
Transport properties show a substantial two- dimensionality of these compounds, due to tt 
states. 

We studied CaSi2 in its trigonal structure, and show that the distortion can be removed 
only by subtracting (in a virtual crystal-like approximation) at least 0.7 electrons. In fact, 
CaBco.sSii.s shows no tendency towards trigonal distortion, and electron states near Ep of 
similar bonding nature to those of MgB2. 
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TABLES 

TABLE I. Phonon frequencies at F for the AlB2type systems (in cm~^). Level degenaracies 
are shown in parenthesis. 





MgB2 


BeB2 


(2) 


328 


339 


(1) 


419 


298 


(2) E2g 


665 


987 


(1) 


679 


639 
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FIGURES 



FIG. 1. Band structure of MgB2. All the energies are referred to the Fermi level, taken as zero. 
FIG. 2. Total and partial density of states (PDOS) for MgB2. 

FIG. 3. Band structure of BeB2. 

FIG. 4. Band structure of MgB2 (panel (a)) and Mgo.5Alo.5B2 at the experimental (panel (b)) 
and frozen MgB2 (panel(c)) lattice parameters. 

FIG. 5. Transport properties of MgB2 and BeB2 as a function of doping in a rigid band scheme. 
Full and dashed lines correspond to BeB2 and MgB2 respectively. Circles (squares) indicate the 
y (z) principal values of the plasma frequency 0,p (in eV), of the mean squared Fermi velocity vp 
(in 10^ cm/sec) and of Rxyz (in 10"^'^ m^/C, notice that Rzxy = Ryzx)- DOS are in states/eV-cell. 

FIG. 6. Band structure of CaSi2 in the ideal AIB2 (dashed lines) and trigonal (full lines) 
structures . 

FIG. 7. Total and partial density of states (PDOS) for CaSi2 in the AIB2 (dashed lines) and 
in the trigonal structure (full lines). 

FIG. 8. Upper pannel: total energy of CaSi2 as a function of the internal parameter z, for two 
different values of pressure, at the corresponding experimental lattice constants ( a = 3.7554 A, 
c = 4.3951 A for P = 15.0 GPa and a = 3.7668 A, c = 4.4752 A for P = 12.8 GPa, after g). 
The experimental values of z are indicated by arrows. Lower pannel: same as above, but with Ca 
substituted by a virtual cation having nuclear charge Z = 20 + x. The larger experimental pressure 
has been used (a = 3.708A and c = 4.028A, the experimantal value of the internal parameter is 
z = 0.448). The zero of energy corresponds to the AIB2 value z = 0.5. 

FIG. 9. Band structure of CaBexSi2_x {x = 1) in the hexagonal Brillouin zone, at experimental 
lattice parameters. 

FIG. 10. Total and partial density of states (PDOS) for CaBexSi2-x {x = 1) in the AIB2 
structures (see text). 
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Fig. 1 G. Satta et al. 





Fig. 3 G. Satta et al. 
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